(N 



A comment on the paper "On the orbit of 
the LARES satellite", by I. Ciufolini 



Lorenzo lorio 

_ Viale Unita di Italia 68, 70125 

, Bari, Italy 

O ; tel./fax 0039 080 5443144 

■ e-mail: lorenzo.iorio@libero.it 



Abstract 



In this note we comment on a recent paper by I. Ciufolini about the 
possibiUty of placing the proposed terrestrial satellite LARES in a low- 
altitude, nearly polar orbit in order to measure the general relativistic 
ly^ , Lense-Thirring effect with its node. Ciufolini claims that, for a depar- 

^ ' ture of 4 deg in the satellite's inclination i from the ideal polar config- 

uration (i = 90 deg), the impact of the errors in the even zonal har- 
monics of the geopotential, modelled with EIGEN-GRACE02S, would 
, be nearly zero allowing for a few-percent measurement of the Lense- 

' Thirring effect. Instead, we find that, with the same Earth gravity 

, model and for the same values of the inclination, the upper bound 

' of the systematic error due to the mismodelling in the even zonals 

O , amounts to 64% of the relativistic effect investigated. 

^ [ Keywords: Lense-Thirring effect; Earth gravity field; polar orbits; new 

satelhte 

; 1 The polar configuration for measuring the Lense- 

5^ ! Thirring effect 

The possibility of measuring the general relativistic gravitomagnetic Lense- 
Thirring effect by means of the node of a LAGEOS-like satelhte placed in a 
relatively low-altitude (a ~ 8000 km), polar {i ~ 90 deg) orbit-POLARES 
in the following— was proposed for the first time by Lucchesi and Paolozzi 
(2001) and, subsequently, criticized by lorio (2002). The benefits of such 
an idea mainly rely in the possibility of using a relatively cheap launcher 
vehicle and in the fact that, for a perfectly polar configuration (i = 90 
deg), all the classical secular precessions induced on the node by the even 
{i = 2,4,6...) zonal (m = 0) harmonic coefficients Ji of the Newtonian 
multipolar expansion of the terrestrial gravitational potential, proportional 
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to cost, vanish. The main drawbacks of such an orbital configuration are as 
follows 

• The satellite's node is perturbed, among other things, by the £ = 
2, m = 1 constituent of the Solar Ki tide whose period is equal to that 
of the spacecraft's node itself: for i 90 deg it precesses very slowly, so 
that Ki would mimic an aliasing secular trend over an observational 
time span of a few years compromising the recovery of the genuine 
relativistic linear trend of interest. This general feature of motion of 
a polar satellite was already recognized by Peterson (1997) and lorio 
(2005a) in the framework of the GP-B mission. 

• This problem is avoided by choosing an inclination a few deg apart 
from the ideal polar configuration. But, in this case, the systematic 
error S/i induced by the mismodelled part of all the even zonal har- 
monics, emphasized by the low altitude of the satellite, is enhanced, 
depending on the accuracy of the gravity model used. lorio (2002) 
used the full covariance matrix of EGM96 (Lemoine et al. 1998) up 
to degree ^ = 20 showing that for orbits just 1 deg apart from z = 90 
deg the impact of the mismodelled even zonals considered amounted 
to 40%. Such an estimate is likely optimistic because of the use of the 
correlations among the solved-for even zonals. 

In (lorio 2002) the possibility of using POLARES in conjunction with the 
existing LAGEOS and LAGEOS II satellites according to the well-known 
linear combination approach was investigated as well: it turned out to be 
unfeasible because of the qTiitc large value of the coefficient with which the 
POLARES node would enter such combinations. 

lorio (2005b) extensively studied the impact of the new Earth gravity 
models by CHAMP and GRACE on the possibility of using a new satellite 
to measure the Lense-Thirring effect. Among other things, the POLARES 
configuration (a = 8000 km and e = 0.04) was re-analyzed with the EIGEN- 
CGOIC model (Reigber et al., 2006), up to degree £ = 20 and, much more 
conservatively than in (lorio 2002), by linearly summing up the absolute 
values of the individual mismodelled classical precessions; the situation is 
now improved with respect to the EGM96 case, but it turned out that for a 
shift of just 2 deg in i with respect to the ideal polar geometry the bias due to 
the mismodelling in all the uncancelled even zonal harmonics still amounts 
to about 25%. In lorio (2005b) also the linear combination scenario with 
LAGEOS and LAGEOS II was investigated showing that for a = 8000 km, 
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e = 0.04, and 60 deg< i <80 deg the systematic error due to the mismodelled 
even zonals is 1-3%. 



2 The departures from the ideal polar configura- 
tion 

The subject seemed, thus, to have exhaustively been treated so far, when a 
new paper on it by Ciufolini (2006) appeared. Basically, the only novelty 
of such work, which mainly reproduces the content of Section 4.2.1 and 
Section 4.4 of lorio (2005b) without quoting it, is a huge underestimation 
of the impact of the uncertainties in our knowledge of the geopotential on 
a certain orbital configuration of POLARES. Indeed, Ciufolini (2006), who 
used the EIGEN-GRACE02S Earth gravity model (Reigber et al., 2005), 
after discussing the problem of the Ki tide, proposed to circumvent it by 
adopting for POLARES an orbital configuration with a = 7878 km and 

1 < 86 deg or z > 94 dog, i.e. an inclination's departure of 4 deg from the 
ideal polar geometry. He explicitly claims that, in this case "[...] if LARES 
would be launched in a nearly polar orbit the use of LAGEOS and LAGEOS 

2 satellites would not be anymore useful in order to reduce the error budget 
(and would indeed only introduce an additional error), since the effect of the 
even zonal harmonics on the node of LARES would be nearly zero, [...]". 
Unfortunately, the situation is quite different, as it could already have been 
inferred from Section 4.2.1 and Figure 5 of (lorio 2005b). For the sake of 
a direct comparison, here we use EIGEN-GRACE02S as well, and in order 
to get a conservative upper bound of the systematic error induced by the 
mismodelling in Ji we linearly sum up the absolute values of the individual 
mismodelled node precessions up to ^ = 40 

40 

Sfi <J2\n.e\sJi (1) 
e=2 

by using the calibrated errors in j£ (Reigber et al. 2005). The coefficients 
Q,£ of the classical node precessions were explicitly worked out up to degree 
^ = 20 in (lorio 2003): for, e.g., £ = 2 we have 

3 I R\ cosi 



where R denotes the Earth's mean equatorial radius and n = -J GMj a? is 
the Keplerian mean motion. Ciufolini (2006) did not explain how he assessed 
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the error due to the even zonals (root-sum-square calculation? Sum of the 
absolute values of the individual errors?), apart from claiming that he used 
the analytical expressions of the nodal precession of a satellite, up to £ = 10, 
from an unspecified reference R. Tauraso (2004). For a = 7878 km, e = 0.04 
and z = 86/94 deg we gelQ 

6fi < 64%. (3) 

In Table [1] we release the details of our calculation. As can be noted, the 
uncancelled precession due to 6J2 amounts to 70% of the entire error. 

In addition to the static part of the geopotential, also its time-dependent 
components must also be considered. In particular, for i = 90 it 4 deg, the 
mismodelled part of the i = 2 m = constituent of the 18.6-year tide would 
have a serious aliasing impact on a sought few-percent measurement, espe- 
cially over an observational time span of just 3 years, as proposed by Ciu- 
folini (2006). The uncancelled secular variations of the even zonals J2,Ja, J& 
would be another source of systematic error. 

Thus, it seems to us very difficult to agree with the conclusion by Ciu- 
fohni (2006) "A nearly polar orbit for LARES at an altitude of about 1500 
km would be suitable for a measurement of the Lense-Thirring effect with 
accuracy of a few percent.". A new satellite can be fruitfully used only in 
conjunction with LAGEOS and LAGEOS II. Such existing satellites, how- 
ever, would set the total realistic accuracy obtainable to a few percent level 
because of the impact of the non-gravitational forces acting on them, inde- 
pendently of how well they could be reduced on LARES. Indeed, it is not 
clear if and how LAGEOS and LAGEOS II could benefit from the reduction 
of the non- gravitational forces on LARES. Such interesting technological 
and engineering efforts (Bellettini et al. 2006) could likely turn out to be 
really and fully useful if the launch of at least two entirely new spacecraft 
was implemented (lorio 2005b). 
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Table 1: Individual mismodelled node precessions 
by the calibrated errors in J^, ^ = 2,4, 6. ..40, in milliarcseconds per year 
(mas yr^^), according to the variance matrix of EIGEN-GRACE02S Earth 
gravity model (Reigber et al. 2005) for a = 7878 km, e = 0.04, i = 86 deg. 
The mismodelled precessions for £ > 30 are smaller than 0.1 mas yr~^. The 
Lense-Thirring effect for such an orbital configuration amounts to 116.6 mas 
yr~^. The upper bound of the total error 5fi < X]^=2 is quoted, in mas 

yr~-^, in the last row: it amounts to 64% of the Lense-Thirring effect. The 
most important contribution comes from J2 whose mismodelled precession 
amounts lo 70% of the total error. 
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0.5 


20 
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22 
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24 
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30 




32 
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36 
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40 




6fi 


68.1 
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